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It has been argued by Zurek and Kibble that the likelihood of producing defects in a continuous 
phase transition depends in a characteristic way on the quench rate. In this paper we discuss 
an improved experiment for measuring the Zurek-Kibble scaling exponent o for the production of 
fluxons in annular symmetric Josephson Tunnel Junctions. We find a ~ 0.5. Further, we report 
accurate measurements of the junction gap voltage temperature dependence which allow for precise 
monitoring of the fast temperature variations during the quench. 

PACS numbers: 11.27.+d, 05.70.Fh, ll.10.Wx, 67.40.Vs 

INTRODUCTION 

The experiments that we describe here are part of a programme 0, 0, designed to see whether continuous phase 
transitions in superconductors proceed as fast as they can, in the sense that the resulting domain structure reflects 
causal horizons. Although, under adiabatic change, correlation lengths do diverge at the critical temperature T Cl in 
reality causality prevents any lengths diverging since transitions take place in a finite time. This seemingly simple 
observation is not of purely academic interest, since causality is universal, and the same is equally true for the early 
universe with its anticipated rich sequence of symmetry breaking. Constraints imposed by causal horizons in the early 
universe are known to have observable consequences Q and, in addition to being of interest in its own right, one 
motive for our, and related, work is to see how transitions in condensed matter systems can mimic this behaviour. 
The argument for the domain structure after a transition being determined by causality was first made by Kibble 
|5j for the early universe and, independently, by Zurek for condensed matter systems, in what we term the 

Kibble-Zurek (KZ) scenario. In particular, if transitions are frustrated, topological defects (vortices, monopoles, etc.) 
arise to mediate the different groundstates of the theory. Their density will be related to the nature of the domain 
structure present and be constrained by causality in turn. The relevance of this is that defects are, in principle, 
readily countable, permitting us to check this proposition. Without taking the parallels any further, we note that all 
reasonable models for the early universe show frustration on cooling 

We would expect that, the faster the quench through the transition, the more defects we would see. Whatever the 
details, and there are several ways to estimate the defect density in the KZ picture, all predict characteristic scaling 
behaviour in the quench time tq (the inverse quench rate) defined by: 
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Specifically, if £ is the separation of defects at the time of their production then it scales with tq as 

where £0 is, most simply, the cold correlation length, and r the relaxation time of the long wavelength modes. 
The scaling exponent a > is, in the mean-field approximation, determined from the static mean-field scaling 
exponents and whether the dynamics is largely wavelike (early universe) or dissipative (condensed matter). That is, 
universality classes of continuous adiabatic transitions lead to identical scaling behaviour of domains as the cooling 
rate is changed. 

For Josephson Tunnel Junctions (JTJs) the topological defect is a fluxon i.e. a supercurrent vortex carrying a single 
quantum of magnetic flux $0 = h/(2e) in the plane of the oxide layer between the two superconductors that make up 
the JTJ. Our experiments, summarised in part in |3(, that we detail below, show that the spontaneous production of 
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fluxons at a temperature quench does, indeed, scale with the quench rate, validating the basic picture of Kibble and 
Zurek. 

This is a considerable achievement since, in general, it has proven surprisingly difficult to establish the scaling be- 
haviour of Eq. (J2J. In fact, although there have been many experiments that have been performed to test the KZ 
picture that are commensurate with it or, when not, explicable in its framework |iA llfxll4llallall4llal 
our experiments are the only experiments with condensed matter systems to date that are sensitive enough to show 
unambiguous scaling behaviour. 

It is worth commenting briefly on how using Josephson Junctions enables us to avoid the main problems that have 
befallen the other experiments with superfluids and superconductors which, superficially, look simpler and more 
direct than ours. The first problem, which besets superfluids, is that of relating the density of defects at the time of 
measurement to the density at the time of formation, as J2J) requires. This was a particular problem for superfluid 
4 He, leading to the null experiments of after the spurious results of Although this is avoided in spontaneous 
vortex production in superfluid 3 He, because of the ability to count vortices more accurately in this case, the fact 
that 3 He is heated by its disintegration under bombardment by soft neutrons 0, ^| means that the fixed rate of 
the nuclear reaction constrains the effective cooling rate. The penalty is that it is not possible to extract scaling 
exponents explicitly. 

Both of these problems look to be avoided wit h p lanar superconductors, where flux is conserved, and where a wide 
range of cooling rates can be implemented |i5j,[l(j. However, since only net flux can be measured, effectively halving 
the scaling exponent, its expected value is so small that it takes an even wider range of quench rates than are available 
in order to show scaling robustly. This is compounded by the fact that, not only are the systems noisy, but Hindmarsh 
and Rajantie 20] have shown that the freezing in of pure magnetic flux gives an additional contribution to the flux 
density anticipated from Eq.Q, which complicates the issue. 

While JTJs possess the virtues of superconductors, they avoid those problems since our samples are so small that 
we expect to see no more than one (conserved) fluxon in each annulus in each measurement. That is, for annuli of 
circumference C, we have C < £. As a result, noise is minimised. Also, because JTJs act as their own thermometers, 
tq can be measured to high accuracy. Further, there is no counterpart of the Hindmarsh- Rajantie mechanism for flux 
produced in a narrow slit. 

However, we do have a potential problem in that Eq.J2J), couched in the language of causal horizons, is not designed for 
systems, such as ours, that are much smaller than them. Instead, we propose that the probability /i for spontaneously 
producing one fluxon in the thermal quench of a symmetric annular JTJ of circumference C < £ should scale with 
the quench time tq (the inverse quench rate) as 



/i-S = ?f^) ■ (3) 



where the scaling exponent a depends on the nature of the junction. We shall provide a justification for J3J in our 
concluding sections but, for the moment, we take it as a consequence of J5J. 

We conclude this section by putting our new experiments in the context of our old. In 2001 our first proof-of-principle 
experiment with JTJs was performed 0, 0| , to test the scaling behaviour of • The experiment consisted of taking 
an annular JTJ isolated from its LHe surroundings and making it undergo a forced phase transition by heating it 
above its superconducting critical temperature and letting it cool passively back towards the LHe temperature with 
no external current or magnetic field. In common with all our experiments, any trapped fluxon can be detected by 
the appearance of a current peak in the I-V characteristic of the JTJ, as will be explained later. 
To derive a two of us had argued earlier [2l| that the relevant causality is provided by the finite velocity of electro- 
magnetic waves in the JTJ, the Swihart velocity [2^,[23j]. Under the idealistic assumptions of a) weak coupling of the 
superconductors and b) exact critical slowing of the Swihart velocity at the critical temperature T = T c , we predicted 
(7 = 0.25|2lj. The experiment was successful, with a commensurate with scaling behaviour © with this value of a. 
The experiments performed subsequently, that we shall describe in this paper, have forced us to revise our assumptions. 
In our recent Physical Review Letters |3| we showed new scaling behavior for the spontaneous production of fluxons, 
in which f\ was indeed seen to clearly follow an allometric dependence on tq, but with a scaling exponent a = 0.5. 
We discuss this experiment in more detail in the subsequent sections and confirm this behaviour with data from new 
samples (see Fig. . The discrepancy between the observations of the early and late experimental values of a may be 
less than it looks at first sight, given the high accuracy of the latter and the relative scatter of the former. However, 
from a theoretical viewpoint, it may also lie in the fact that we need to take into account the consequences of finite 
size and the nature of the fabrication of the junction, which differs between the two sets of samples. In this regard 
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realistic condensed matter systems cannot match the early universe for their extension and uniformity, upon which 
P)l. with a based simply on the usual critical exponents, is predicated. Again we postpone a proper discussion of 
these issues, and the calculation of a, to the concluding sections of this paper. 

However, before then we shall describe the nature of the samples and the experimental setup, the use of the JTJ as 
its own thermometer to measure quench times, and then the results confirming the scaling behaviour of @. 




FIG. 1: Layout of of the 4.2mm x 3mm Si chip containing four series biased Nb — Al / Al ox / Nb Josephson tunnel junctions. 
It integrates three ring shaped junctions having a mean circumference C = 500 fim and a width Ar = 4 fim, one 4 x 500 fim 2 
overlap-type linear junction and two meanderline resistive Mo strips used for heating. 



THE SAMPLES 

To begin with some generalities, the annular JTJs (AJTJs) that we have used are high quality Nb/Al — Al ox /Nb 
JTJs fabricated on silicon substrates using the trilayer technique in which the junction is realized in the window 
opened in a SiO insulator layer. The so called "idle region", i.e. the overlapping of the wiring layer onto the base 
electrode was about 3 \im for all the junctions. The thicknesses of the base, top and wiring layer were 200, 80 and 
400 rim, respectively. Details of the fabrication process can be found in Ref.[36j. The samples were fabricated at the 
Superconducting Electronics Laboratory of the Institute of Radioenginecring & Electronics of the Russian Academy 
of Science in Moscow, while the measurements were carried out at the Physics Department of the Danish Technical 
University in Lyngby(DK). All the experiments performed to date have been carried out on AJTJs with a mean 
circumference C — 500 fim and a width Ar = 4 fim. AJTJs with larger circumferences have been fabricated, but we 
have yet to use them. 

The new chip layout developed for the K-Z experiment is shown in Fig.^ It integrates three ring shaped junctions 
having a mean circumference C — 500 fim and a width Ar = 4 \xm and one 4 x 500 fim 2 overlap-type linear junction. 
The four JTJs are biased in series. The rightmost A JTJ was obtained by the superposition of two superconducting 
rings, as depicted in Fig. while the two AJTJs in the layout middle were realized by the superposition of a ring 
shaped top electrode over a superconducting plane, as shown in Fig. ^jp. Here we anticipate that this difference in 
the sample topology did not affect the measured spontaneous defect production. 

For the new experiment a faster and more reliable heating system was required. This was achieved by integrating a 
meander line 50/im wide, 200nm thick, and 8.3mm long Mo resistive film in either ends of the 4.2mmx3mmx 0.35mm 
Si chip containing the Nb/AlOx/Nb trilayer JTJs (also shown in Fig. QJ. These resistive elements have a nominal 
d.c. resistance of 50f2 at LHe temperatures and, due to their good adhesion with the substrate, are very effective in 
dissipating heat in the chip. In fact, voltage pulses a few /is long and a few volts high applied across the integrated 
heater provided quench times as low as 1ms, that is more than two orders of magnitude smaller than for the previous 
situation The on-chip heaters are able to sustain thousands and thousands of pulses as well IV continuous bias 

without any appreciable change of their electrical resistance. 
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(a) (b) 

FIG. 2: Sketch not in scale of two geometries for the annular Josephson tunnel junctions used for the Kibble-Zurek experiment. 
The junction base and top electrodes a shown in dark and light gray, respectively. In (a) the annular junction is obtained by 
the superposition of two superconducting rings, while in (b) it is realized by the superposition of a ring shaped superconducting 
top electrode over a superconducting plane base electrode. 



THE EXPERIMENTAL SETUP 

A detailed description of the experimental setup has been given in Rcf . 2] . Briefly, the chip with the A JT Js is mounted 
to a Cu block enclosed in a vacuum-tight can immersed in the liquid He bath. A Ge thermometer anchored to the Cu 
block allowed for precise measurements of the chip static or slowly changing temperature. The junction itself is used 
for quickly changing temperature measurements, as it will be explained in the next section, but with lower accuracy. 
The chip was heated above the AJTJ critical temperature by a voltage pulse applied to one or both of the integrated 
meanderline heaters. Then the heat is removed from the chip both through the thermal contact with the Cu block 
and by He exchange gas inside the can. For all the experiments performed so far the He exchange gas pressure was 
fixed at a value of about 7mbar which makes the heat flow through the copper base plate predominant. 
In order to minimize thermal gradients during the thermal cycles, particular care was taken to have clean and flat 
contact surfaces of both the chip and the Cu block; further, the voltage pulses were applied simultaneously to both 
the on-chip heaters. A tiny layer of low temperature grease was span underneath the chip to improve the thermal 
contact. 

At the end of each cycle the possible spontaneously generated fluxons are static. An external current supplied to 
the AJTJ sets the fluxons (if any) in motion around the annulus and quantized voltages develop across the junction 
itself. In other words, we count the number of produced fluxons by a careful inspection of the junction current-voltage 
characteristic(IVC). (Due to the annihilation of a fluxon-antifluxon pair, this idea works well as long as the chances 
to spontaneously generate two fluxons are small.) Fluxon motion at 4.2K is very unstable in our samples due to very 
low junction losses, therefore the IVC was better inspected at higher temperatures where larger losses stabilize the 
fluxon motion. For our sample the optimal temperature to look for fluxons moving around the ring was in a range 
between 6 and 7K. 

Quenching experiments were carried out in a double /x-metal shielded cryostat and the sampleholder can was sur- 
rounded by a superconducting Pb shield. In turn the chip holder inside the can was enclosed in a in a cryoperm 
shield and in one more superconducting Pb shield. The transitions from the normal to the superconducting states 
were performed with no current flowing in the heaters and the thermometers. During the quench the JTJ was also 
electrically isolated: in fact, both the junction voltage and current leads were open during the whole thermal cycle. 
Finally, all the measurements have been carried out in an electromagnctically shielded environment. 

In order to run batches of several thousand equal thermal cycles with given parameters, automatization of thermal 
cycles was implemented by means of a switching unit controlled by a GPIB interface; that also allowed for much 
more robust statistics to be achieved. At the end of each thermal quench the junction IVC is automatically digitally 
acquired, converted and stored. Then at the end of each batch with a given value of tq an algorithm has been developed 
for the analysis of the large amount of IVCs and the automatic detection and counting of the spontaneously trapped 
fluxons. 



5 



DETERMINING THE QUENCHING TIME 



The quench time tq was continuously varied over more that four orders of magnitude (from 20s down to 1ms) by 
varying the width and the amplitude of the voltage pulse across the integrated resistive elements. In order to estimate 
the quenching time tq we use the observation by Thouless |34j that the junction itself acts as a thermometer, as far 
as its temperature is below the critical temperature. This gives us an unrivalled accuracy in measuring tq over other 
experiments looking for KZ scaling behaviour. Specifically, the temperature dependence of the gap energy A in a 
strong-coupling superconductor, such as Nb, 



ACT) , A(T)T C 

also applies to the junction gap voltage V g that is proportional to A. An exp erimental demonstration of Eq.Q in 
Nb/Nb tunnel junctions having native oxide was first evidenced in 1976 |35| . well before the development of the 
trilayer technique used here that, by exploiting the more compact Al oxide, allows for a higher quality and a more 
stable tunnel barrier. 

Eq.|@J can be manipulated: considering that arctanhx = \ In as far as x 2 < 1 and approximating ln(l ± x) ~ 
±x — x 2 /2 ± £ 3 /6, we easily get: 



A(T) 
A(0) 



-i 2/3 



In our samples we have found that, provided the JTJ is current biased at about 20-25% of the total current jump at the 
gap voltage, the junction temperature could be monitored to a high degree of accuracy and speed by resorting to Eq.Q). 
This is shown in Fig. |21 where the experimental values (open circles) of the junction gap voltage V g (T) = 2A(T)/e 
are plotted at different values of the temperature as measured by a calibrated Ge thermometer. The solid line is the 
prediction that follows from Eq.Q with fitting parameters T c = (9.12 ± QM)K and V g (0) = (2.89 ± 0.02)mF, i.e. 
2A(0)/fcsT c ~ 3.73 ± 0.03. As a result, Eq.QJ can be used efficiently to estimate the junction temperature to a high 
degree of accuracy for V g (T) > lmV, i.e., T < 8.5K. Above 8.5K the experimental data for the gap voltage saturate 
to a finite temperature-independent value corresponding to the product of the JTJ normal resistance and the bias 
current. This way of monitoring the system temperature is particularly convenient when the temperature changes 
rapidly and a complete thermal cycle occurs on a ms time scale or even shorter. However, the overall temperature 
accuracy ST = \-^r\5V g cannot be smaller than 2mK due to a voltage accuracy 5V g of about 2[iV on a fast digital 
oscilloscope. 

Assuming the chip exchanges heat mainly through a massive copper base plate (with thermal constant t±) which 
in turn exchanges heat with the surrounding helium gas (with thermal constant t 2 ), the thermal relaxation during 
the thermal quench has been fitted by a double exponential decay of the form: 



T(t) = T fm + AT, exp f-^p) + AT 2 exp (-^^) > (5) 

with ATi, AT 2 , T\ and r 2 fitting parameters. Tf in and to are known from the experiments. If the time origin is 
triggered by the voltage pulse, then t corresponds to the time at which the pulse ends. Once the parameters in 
Eq. (JjJJ are determined or measured, the quenching time tq can be inferred from its definition Eq. JTJ , after Eq. JSJ) has 
been extrapolated up to the critical temperature Tc- At the end of this process of fitting and extrapolation, tq is 
known to an overall accuracy of about ±10%. 



THE K-Z MEASUREMENTS 



The experimental results reported here refer to three identical AJTJs belonging to two different chips made within 
the same batch having a critical current density J c (0) ~ 60A/cm 2 yielding a Josephson penetration depth Aj(0) ~ 
50/im. Assuming a' = a = 3.5 in Ea. (|ll|) . this leads to a value of £0 ~ 17/im. For all samples the high quality has 
been inferred by a measure of the I-V characteristic at T = 4.2 K. Due to the very high reliability of the fabrication 
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FIG. 3: The temperature dependence of gap voltage V g . The open circles are the experimental data with the junction biased 
at 25% of the total current jump at the gap voltage; the solid line is the prediction that follows from Eq.JIJ. The fit yields 
T c = (9.12 ± 0.04)K and V g (T = 0) = (2.89 ± 0.02)mV. 



line the sample tunnel barriers have the same geometrical and electrical parameters. However, to distinguish them we 
will name them after their chip and junction numbers, that is, we have samples 08-3, 08-4 and 22-4. Samples 08-3 and 
08-4 belonged to the same chip. Samples 08-4 and 22-4 had the geometry sketched in Fig.|2K, while sample 08-03 had 
the geometry shown in Fig. 12). The solid symbols in Fig.0]show on a log-log plot the measured frequency fi = n\/N 
of single fluxon trapping, obtained by quenching the samples N times for each value of a given quenching time tq, ri\ 
being the number of times that the inspection of the low temperature AJTJ current- voltage characteristics at the very 
end of each thermal cycle showed that one defect was spontaneously produced. N ranged between 100 and 5000 and 
n\ was never smaller then 10, except for some of the rightmost points (tq > 10s)for which ri\ > 3. All samples has 
undergone a total of more than 100,000 thermal cycles without any measurable change of their electrical parameters. 
The vertical error bars gives the statistical error j\j\Jn\. The relative error bars in tq amounting to ±10% are as 
large as the dot's width. 

It is quite evident that the dependence of the trapping frequency on the quenching time is the same for all three 
AJTJs independently on the geometry of their base electrodes. Careful measurements of the junction IVC during the 
N-S transition indicate that the junction critical temperature might differ from the critical temperature of the base 
electrode film by no more than lOmK, while the critical temperature of the wiring film, being twice thicker, exceeded 
the critical temperature of the base electrode by about lOOmK. In other words, at the time the Josephson effect is 
installed, the base electrode is only weakly superconducting and cannot exercise any shielding effect. This explains 
why we have not observed any difference between the data of junction 08-3 and those of junctions 08-4 and 22-4. 
To test Eq.(|3Jl, we have fitted the data of all samples with the same allometric function fi = otq 6 , with a and b as 
free fitting parameters. A linear regression of log/i vs. log tq, represented by the continuous line in Fig. 01 yields 
a = 0.01 ± 10% (taking tq in seconds)and b = 0.51 ± 5%. Therefore the present experiment suggests that the scaling 
exponent is a = 0.5, rather than the value a = 0.25 suggested by our earlier attempt^,)^. For comparison, the data 
of this experiment are reported in Fig. ^ as open stars. 

The shift in intercept (or, equivalently, prefactor a) between the two sets of data is to be expected. The AJTJs 
used in the first experiment, although of the same geometry as samples 08-4 and 22-4 (see Fig. Et), had a Josephson 
current density J c about 60 times larger. As we shall reiterate later, this means a smaller Zurek length £ oc l/y/j c , 
with a correspondingly greater likelihood of observing a fluxon. In other words, we have moved from a situation in 
which C < £ to C <C £. As we shall see, the value a — 0.01 ± 10% (taking tq in seconds) is 6-7 times larger than the 
predicted value. As a bound we only expect agreement in the overall normalization a to somewhat better than an order 
of magnitude. Empirically, the different condensed matter experiments have shown that the ratio a bserved/&pr£djfted 
varies widely from system to system; 0(1) for superfluid 3 He 0,0], very small for high-T c superconductors |l6| . 
Although the best fit to the stars alone is a — 0.25, when seen in conjunction with our new data, this value is not 
so compelling, due to the poor statistics and to the scattering of the stars. In retrospect, we do not exclude the 
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possibility of systematic error in the 2001 experiments. The trapping frequency /i increases when the thermal cycling 
occurs in an externally applied magnetic field. Therefore an insufficient shielding of the earth's magnetic field might 
be the cause of the systematic shift upwards of /i for increasing tq in the first experiment. For the actual experiment 
we have taken a lot of precautions against the possibility of having a significant d.c. residual magnetic field. In fact, 
although its absolute value cannot be measured, we have checked that the measured trapping frequency did not change 
i) by rotating the sampleholder in the horizontal plane (which changes the sample orientation with respect to the 
direction of the earth's magnetic field) and ii) by rotating the chipholder inside the sampleholder, which is kept fixed 
(which changes the sample orientation with respect to the shields). Furthermore, detailed measurements have been 
carried out of the dependence of the trapping frequency on the strength of an external field applied perpendicular to 
the junction plane for different samples and for several values of the quenching time. Such data will be reported in 
another paper |3^| that will also discuss the interplay between the Abrikosov vortices produced in the superconducting 
junction electrodes and the Josephson vortices produced in the junction barrier during the thermal quench. As far as 
the measurements presented in this work are concerned we made sure that for each value of the quenching time the 
corresponding trapping frequency lies at the minimum. 

There is another plausible explanation for the different experimental findings to which we shall return in the next 
section. This is that the samples used in the two experiments, although constructed in the same fabrication line, 
might have slightly different proximity interfaces, which results in quite different temperature dependence of the 
critical current density near the transition temperature. 

Whatever the reasons for the differences, the new data of Fig. 01 resolve another issue. It could be argued that what 
we are really seeing is phase ordering from the individual superconductors, in which case /i would only be a measure 
of the flux trapped in just one of the two superconducting annuli forming the AJTJ, not including those cases in 
which flux is trapped simultaneously in both rings. If we had taken the relevant velocity to be that of phase ordering 
in the individual superconductors, as invoked by Zurek 0, Q when considering the spontaneous flux produced on 
quenching annuli of simple superconductors, at the same level of approximation we would have predicted a = 0.25, 
which is manifestly not the case. 

There is, however, one unexplained observation. Unlike in our earlier experiment, in this experiment we never 
observed the simultaneous production of two (or more) fluxons (or antifluxons), despite the fact that we have detected 
more than 10, 000 single fluxons for each sample. This is in contradiction with the expectation that, extrapolating 
from an infinite junction, the probability fa of trapping two homopolar fluxons is f% = j\ /2, as far as the productions 
of a defect in different ring regions are considered independent events. 




FIG. 4: Log- log plot of the measured frequency /i of trapping single fluxons versus the quenching time tq. Each point 
corresponds to many thermal cycles, closed squares for sample 22-4, closed triangles for sample 08-3 and closed circles for 
sample 08-4. The vertical error bars gives the statistical error while the relative error bars in tq amounting to ±10% are as 
large as the dots' width. The solid line is the fit of all data to an allometric relationship fi = a 6 which yields a = 0.01 ± 10% 
(taking tq in seconds) and b = 0.51 ± 5%. The open stars represent the data obtained in a previous experiment 0,0- 
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THEORY 

The value of a — 0.5 is obviously in disagreement with our earlier prediction of a = 0.25 given in |21| . Two possible 
reasons need to be explored. The first is that our assumptions of idealised JTJs with weak coupling and exact critical 
slowing down of the Swihart velocity at the transition used in [2l| are not valid. The second is that the use of causality 
seems suspect when dealing with systems like the annuli discussed here that are smaller than the KZ causal length £. 
That is, how should we treat causal horizons that are larger than the systems at the time that defects are formed, if 
we are using them to define defect separation? Is it the case that Eq.© follows from Eq.Q? 

We address the second problem first. In seeing how Eq.(jHJ can arise we are helped by the fact that, since the 
original KZ papers were published, there has been considerable analytic and numerical work performed for ideal 
systems obeying dissipative equations EE EE EE EE EE EE EE EE| , m particular the time-dependent Ginzburg- 
Landau equation. Although there has been no attempt to model Josephson tunnel junctions we can draw several 
conclusions for JTJs from the phase transitions in simple systems, such as superfluids and superconductors, that 
suggests that small system size is not aproblem in principle. This suggests a rather different picture from that 
originally posed by Kibble and Zurek in [EIEQ]) but with the same outcome in (01 for large systems. 

As originally posed, the separation of defects is directly related to the correlation length. By definition, the 
correlation length is determined from the large distance behaviour of the correlation function (and thereby the position 
of the nearest momentum-space singularity of the power spectrum in the complex plane) . On the other hand, analytic 
and numerical work shows that defect separation is more strongly related to the short-distance behaviour of the 
correlation function, since it is this that controls the field zeroes (in this case modulo 27r) with which simple defects 
are associated j.Slj . In the first instance after the transition this is given by the moments of the power spectrum, rather 
than its singularities. Of course, if there is only a single length in the model at the time the transition is implemented, 
the two must be essentially identical, although they play very different roles [30|. Provided thermal fluctuations 
are controlled, this is plausible for a dissipative system [32j. Field zeroes, which can mature into defects, initially 
occur on all, particularly small, scales. Most annihilate very rapidly. What drives those that have not annihilated to 
become the cores of defects on large scales is the unstable growth of long wavelength modes, which transfers power to 
long distance fluctuations at the expense of small. With fluctuations starting from such small beginnings, from this 
viewpoint the fact that, ultimately, £ > C, does not hinder defect formation. 

Suppose now that we increase C so that C > £ and we see fluxons every time. f\ is then not a useful measure and, 
instead, we measure total flux i.e. the variance An in the net number n of fluxons (i.e the number of fluxons minus 
the number of antifluxons) . Using the spacing of Eq. (J2J a random walk in phase along the annulus suggests 

\2 _ c _ C ( T Q 



Since (n) — 0, (An) 2 = (n 2 ) = /i when f% < 1 is sufficiently small that we can neglect the possibility of seeing 
more than fluxon. As a result we rederive Eq.ljJJ, but without having had to invoke causality in the same way. More 
work is being done, particularly on numerical simulations with periodic boundary conditions |4oj |. but for the sake of 
argument we assume that finite size is not the reason for the disparities in a. 

It is a second issue as to whether the time scale for defect formation, the time it takes for long wavelength modes 
to evolve fully, is related to the KZ causal time at which defects are formed, from which J5J follows. However, very 
simple arguments 0> EE EE show that this is the case, for relatively weak couplings at least, up to logarithmic 
corrections. This enables us to use the KZ scenario to calculate a even though causality is not, directly, the driving 
mechanism. 

We now turn our attention to the specific properties of the JTJs. To reiterate, the theory in Ref.[2l| was developed 
for JTJs whose electrodes are weak coupling superconductors for which the temperature dependence of the critical 
current density J C (T) is given by the Ambegaokar-Baratoff eauation|37): 

^(T) = ^tanh^, (7) 
2 epN 2kbT 

where A(T) is the superconducting gap energy and Pat is JTJ normal resistance per unit area. Eq.Q provides a 
linear decrease of J c near T c : 

Jc(T)=aJ e (0)(l-P), (8) 
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in which the dimensionless quantity a is approximately equal to 2A(0)/fcsT(7 = 3.5. However, our JTJs are based 
on Nb, a strong- coupling superconductor, for which Eq.JTJ) is not necessarily valid. In practice, high quality and 
reproducible barriers are achieved by depositing a thin Al overlay onto the Nb base electrode which will be only 
partially oxidized, leaving a Nb — Al bilayer underneath having a non-BCS temperature dependence of the energy 
gap and of the density of states . The proximity effect in Nb — Al/Al ox /Nb JTJs has been extensively studied and 
it is known to influence the electrical properties of the junctions, such as the current-voltage characteristic and the 
temperature dependence of the critical current density. Specifically, the proximity effect in S — N — I — S junctions 
is responsible for an otherwise subdominant temperature dependence of the critical current density |ssj dominating in 
the vicinity of T c as: 



J c (T)~a'J c (0)(l-;|) , (9) 

where a' is a constant depending on the degree of proximity. The last equation models the tail shaped dependence of 
J c vs. T near T c ; it has been theoretically derived and experimentally confirmed by Golubov et al. 39] in 1995. 
Assuming that the proximity effect is important for our samples here then, on rephrasing the arguments of Ref. |2lj 
with Eq. © replacing Eq. JHJ) , the Josephson penetration depth A j (T) , which plays the role of the system equilibrium 
coherence length £(T), now diverges linearly near T c as: 



e(T) = Xj(T) c< J c (T)-^= Co (l - |) 



(10) 



where 



iQ y 2enod s a' J c (0) ' 

d s being the electrode thickness . 

At time t close to the transition, the temperature T(t) satisfies 

i-saUA (12 , 

1 c J TQ 

where T(0) = T c . The first assumption in the KZ scenario is that causality establishes a time t at which domains and 
defects form, at which time defect separation is £(T(t)). For JTJs fluxon separation is, from l)l()[l. 

Z = Z(T(i)) = \j(T(t))=Z ?jl, (13) 

instead of the £0 y/Tq/t behaviour of [2lJ that we used in Q, ■ 

We have argued that, although it is instability, rather than direct causality, which drives scaling behaviour, it arises 
in a way that is quantitatively indistinguishable from the KZ scenario. 

Thus, even where causality is an inappropriate mechanism we can still adopt its results, that the earliest possible 
time t at which defects could possibly appear satisfies 



at) = -g(i), (w) 

where t is the causal time and c is the Swihart velocity. 

As we said, in Ref.|2l| we had assumed that the Swihart velocity vanished at T c whereas, for realistic JTJs, it just 
becomes very small. Swihart has demonstrated that for a thin-film superconducting strip transmission line the 
solution for the velocity varies continuously as one passes through the critical temperature into the normal state. As 
a result, we assume c(t) — c nn near the transition temperature where c nn is the speed of light in a microstrip line 
made of normal metals. In the case of a microstrip line made by two electrodes having the same thickness d s and the 
same skin depth (5, with d s << <5, separated by a dielectric layer of thickness d ox and dielectric constant e, its value, 
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2 / d ox d s 



depends on the temperature very weakly, but depends on the frequency / through 5 — y p/irfif, p being the normal 
metal residual resistivity. 

The match c(T c ) = c nn is certainly realistic and we still have approximate critical slowing down insofar as c nn is 
much smaller than the zero temperature Swihart velocity Co = \J d ox /2\Loe/i, i.e., when the zero temperature London 
penetration depth A^o << S 2 /d s . For 300nm thick Nb electrodes (p = 3.8^f2cm and Alo = 90nm), 8 ~ 1mm at 
say / = 10kHz, so the last inequality is fully satisfied. At the same frequency, for a value of the specific barrier 
capacitance c s = e/d ox = 0.02F/to 2 typical of low current density Nb — Al/Al ox /Nb JTJs, we get c nn = 7 ■ 10 3 m/s 
and Co = 1.4- 10 7 m/s. We stress that c nn is finite and order of magnitudes smaller than the zero temperature Swihart 
velocity cq. 

The solution of the causality equation Eq. 1|14|) with a non- vanishing Swihart velocity yields a new expression for the 
so called Zurek or freeze out time t — \/ £aTQjc nn = \f T o T Q, where To = £o/c~, m (to = 0{lns)). 
Inserting the value of t into Ea. H10(l we obtain the new Zurek length £: 

i=m = \/^rQC n n=^[S) ■ (15) 

We reach the important conclusion for realistic JTJs that the probability f\ for spontaneously producing one fluxon 
in the quench is still predicted to scale with the quench time tq according to Eq. J3J , but the critical exponent is now 
a = 0.5, rather than a — 0.25. 

By varying tq in the experimentally achieved four decade range 1ms < tq < 10s, we get 10/is < i < Ims that 
is always much larger than To; it means that by the time the Josephson phase 'freezes' the Josephson effect is well 
established. 

The freezing temperature T = T{t) — T c { \ — t / tq) comes out to differ by the critical temperature itself by an amount 
T c — T = T c ^/tq/tq. Further, in the same tq range the normalized freezing temperature T/Tc at which the defect is 
formed is 0.99 < T /Tc < 0.9999. It would be really hard, if not impossible, to measure the temperature dependence 
of J c and c so close to Tc- We have then to resort to theoretical predictions, that is to Eq.@ for the temperature 
dependance of the Josephson current density J C ()T. Further, according to Swihart's calculations and figures |22j. we 
are allowed to assume c(t) = c nn near the transition temperature. 



CONCLUSIONS 



Eq.© is amenable to further experimental tests with A JTJs having different critical current densities J c (0) and/or 
circumferences C. Such experiments should still show the critical exponent a = 1/2, however the prefactor should 
change accordingly. A test of the importance is provided by working with asymmetric Nb — AIN — NbN junctions. 

According to the theory presented in this paper, the expected dependence of the single fluxon trapping frequency 
on the critical current is weak, being proportional to its fourth root, but critical current densities 20-30 times larger 
should produce detectable effects. The change in trapping frequency will be more easily to observe with larger diameter 
annuli due to the linear dependence on C in Eq.©. 

It is worth commenting here on the effect of the unavoidable thermal gradients in physical system undergoing a 
thermal quench; the Zurek-Kibble scaling law Eq.© was derived assuming that thermal gradients are not a problem. 
More precisely, according to the general theory of Ref.^3|, the maximum thermal gradient VT* allowed across a given 
system undergoing a thermal quench at the time of the transition amounts to (T — Tc)/£. Below this threshold a 
saturation is expected in the spontaneous production of defects. In particular, in our case, where tq ranges in the 
interval 10s — 1ms, we have VT* = Tc/(c nn TQ) ranges from 6.5pK/mm to 65mK/mm across our 160/xm diameter 
ring corresponding to a critical value for the maximum temperature difference AT max — VT* x C/ir i=s lpK — lOmK. 
As stated before, our voltage accuracy in the measurements of the gap voltage in contiguous JTJs on the same 
chip allows us to resolve temperature differences as small as few millikelvin; the simultaneous measure of the time 
dependence of the gap voltages of contiguous junctions in response to the shortest heat pulse (i.e. corresponding to 
tq = 1ms) revealed that the temperature difference, if any, is less than lOm-fT. Further, the fact that we have not 
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observed any saturation of the defect production in a wide tq range should indicate that thermal gradient still do 
not affect the defect production in our experiments. Again, in future experiments devised to understand the role of 
the thermal gradients, they can be enhanced by progressively reducing the thermal coupling between the chip and its 
surroundings and/or by supplying the voltages pulses only to one of the two resistive stripes integrated at the chip 
extremities. On the opposite, if needed, it is also possible to further reduce the thermal gradients by using an annular 
integrated heater superimposed on and concentric to the annular JTJ. 

In summary, we see this experiment as providing unambiguous corroboration of Kibble-Zurck scaling over a wide 
range of quenching time to in accord with our predictions for Nb — Al/Al ox /Nb JTJs. As such, it replaces the 
experiment reported in 0, Q by being more realistic theoretically and more sophisticated experimentally. 

We stress that to date, we are the only group to have performed experiments on JTJs and this experiment is the 
only one to have confirmed the Zurek-Kibble scaling exponent for a condensed matter system. We note, however, 
that scaling has been observed 01 m non-linear optical systems, which satisfy equations of the time-dependent 
Ginzberg-Landau type. 
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